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TITLE: FAST COMPUTATION OF DECISION FEEDBACK EQUALIZER 
COEFFICIENTS 

SPECIFICATION 
CROSS REFERENCE TO RELATED APPLICATION 

This application claims priority to U.S. Provisional Application Serial No. 60/322,994, 
filed September 18, 2001, which is incorporated herein by referoice in its entirety. 

1. FIELD OF THE INVENTION 

This invention relates generally to digital communications; and more particularly to 
decision feedback based equalizers. 

2. BACKGROUND OF THE INVENTION 

The structure and operation of communication systems is generally known. Many 
communication systems carry data, e.g., voice, audio, video, file, or other digital data that is sent 
from a transmitter to a receiver. On the transmitter side, data is first formed into packets. This 
data may be raw data or encoded data that represents the raw data. Each of these packets also 
typically includes a header, a known training sequence, and a tail. These packets are then 
modulated into symbols and the symbols are transmitted by tiie receiver and intended for receipt 
by the receiver. The receiver then receives the symbols and attempt to extract the data from the 
packets that are carried by the symbols. 

A "channel" carries the symbols firom the transmitter to the receiver. The channel is 
serviced by a wired, wireless, optical, or another media, depending upon tiie communication 
system type. In many communication systems, such as terrestrial based wireless communication 
systems, satellite based communication systems, cable based communication systems, etc., the 
channel distorts tiie transmitted symbols, fix)m flie perspective of tiie receiver, causing 
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interference between a subject symbol and a plurality of symbols surrounding the subject symbol. 
This type of distortion is referred to as "inter-symbol-interference" and is, generally speaking, the 
time-dispersed receipt of multiple copies the symbols caused by multipath. The channel also 
introduces noise into the symbols prior to their receipt. Each of these concepts is well known. 

Equalizers are now generally employed in an attempt to remove channel effects ftom a 
received symbol stream. Thus, equalizers are essential building blocks of modem receivers, 
especially in broadband applications where inter-symbol-interference is a critical problem. In a 
typical equalizer, the channel between the transmitter and the receiver is first estimated based 
upon the training sequence contained in one or more preambles. Then optimal equalizer 
coefficients (also referred to as taps and/or tap coefficients for the equalizer) are estimated based 
upon the channel estimate. The optimal equalizer coefficients are then used by the equalizer in 
extracting the data from the packet. The optimal equalizer coefficients may also be computed 
after extraction of the data from the equalized data stream based upon blind channel estimates. 
Equalizer coefficient generation should be repeated as often as possible, especially in fast varying 
channel cases, to generate new equalizer coefficients. The received data stream is usually 
buffered during the period that is required for channel estimation and equalizer coefficient 
computations. Thus, the preamble (and also actual data) contained in a packet may be used to 
generate the channel estimate and optimal equalizer coefficients that are employed by the 
equalizer to extract the data from the packet. 

As symbol rates increase and modulation schemes become more complex, equalizers 
have increasingly greater importance. A critical factor in increasing the effectiveness of these 
equalizers is tiie complexity of optimal equalizer coefficient computation. A reduction in this 
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complexity: (1) reduces the memory size required to buffer the received symbol stream sequence 
during the period required for coefficient computations; (2) allows more frequent uploading of 
new coefficients thus enabling the equalizer to track fast channel variations; and (3) simplifies 
the hardware and, resultantly, the die area required for coefficient computation. 

FIG. 1 is a block diagram illustrating a discrete time symbol-spaced Decision Feedback 
Equalizer (DFE) based channel equalization model 100. The channel equalization model 100 
includes a channel 102, a Feed Forward Equalizer (FFE) 104, a Decision block 106, and a Feed 
Back Equalizer (FBE) 108. An input sequence x(n)is complex, independent and identically 
distributed with unit power. Additive noise v(«) is white Gaussian with power spectral density 
. Furthermore, the decisions x(n - d) are assumed to be correct, and hence equal to x(n - S) . 
This assumption makes the design of the FBE 108 and FFE 104 easier, but at the expense of 
introducing error propagation due to possibly wrong decisions. The FFE 104 fimction G(z) has 
length L. The channel (impulse) response vector of the channel h is given in Equation (1) as: 

h = [h{0) h(\) ••• hiN-l)] Equation (1) 

The number of coefficients (taps) M of the FBE 108 function B(z) is assumed greater or equal to 
the chamel memory, i.e., M>N-l. These modeling assumptions are feasible in practice. 

M estimating FFE 104 and FBE 108 equalizer coefficients, the goal is to minimize the 
mean square error quantity of Equation (2). 

C = ^x{n-S)-x{n-^\ Equation (2) 

where x{n-S)is the delayed input signal estimate prior to the Decision block 106. By 
collecting the coefficients of both G(z) and B(z) into vectors, we can express the received signal 
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jc(« - S) in Equation (3) as: 



x„ = y„g-x„b 



Equation (3) 



A channel output model defining y„ may be expressed by: 

y„=x„H+v„ Equation (4) 

where ^ is the (N+L-l)xL convolution matiix corresponding to the channel response and 
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Equation (5) 



In this model, x„ is the 1 x (iV + Z - 1) input vector. 



x„ = [x{n) x{n-\) ••• x(n-N-L + l)] 



Equation (6) 



j;„ is the Ix L input regression vector to the FFE 104, 



y, = [yiri) y{n-\) - x-yin-M)] 
X is the 1 X M input regression vector to the (strictly causal) FBE 1 08, 



Equation (7) 



x„ = [x(n-d-l) x(n-S-2) ••• x(n-S-M)] 



Equation (8) 



and v„ is the 1 x L vector noise process. 

The current effidait methods for computing the optimal filter coefficients of a decision 
feedback equalizer, which optunizes (1), are based on tiie well-known Cholesky decomposition 
method (fi-om a finite-dimension problem formulation). Two published papers: (1) N. Al- 
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Dhahir and J.M. Cioffi, "MMSE Decision-Feedback Equalizers: Finite-Length Results," IEEE 
Trans, on Information Theory, vol. 41, no. 4, pp. 961-973, July 1995; and (2) N. Al-Dhahir and 
J.M. Cioffi, "Fast Computation of Channel-Estimate Based Equalizers in Packet Data 
Transmission," IEEE Trans, on Signal Processing, vol. 43, no. 11, pp. 2462-2473, Nov. 1995 
provide one procedure for computing optimal DFE settings. These equations are referred to 
hereinafter as "Al-Dhahir's equations." 

Generally speaking, Al-Dhahir's equations, as well as oiher existent techniques rely on the 
use of the generalized Schur algorithm for fast Cholesky decomposition of the matrices mvolved 
in both the FBE and FFE optimal coefficient setting computations. However, the overall 
procedures for calculation of the DFE (FBE and FFE) coefficients have the following problems: 

1. These procedures require the use ofnonstructured recursive equations. These equations 
are difficult to implement from the perspective of integrated circuit design. In particular, 
the recursive equations used in Al-Dhahir's DFE tap (coefficient) computer often requires 
the use of a DSP processor. As is generally known, the use of a DSP processor in a real- 
time communication system application is severely constricts system throi^put. 

2. These procedures, in particular Al-Dhahir's equations require a complex tap computer. 
In other words, they require the use of a relatively large number of complex multiplies. 

3. These prior art DFE coefficient computation techniques can be used only for equalizers 
that use a fixed equalizer delay {d), which is set to its maximum value L-l, where L is 
the length of the FFE 104. The prior art techniques cannot use a different delay S of the 
equalizer. 

Thus, there is a need in the art for a new equalizer coefficient computation methodology 
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and system in which such methodology may be implemented. 

SUMMARY OF THE INVENTION 

In order to overcome the above-cited shortcomings, among otiiers of the prior art, the 
method and system of the present invention computes optimal Decision Feedback Equalizer 
(DFE) coefficients J from a channel estimate h. As contrasted to the prior art, the 
methodology of the present invention relies on iterating simple structured recursions to yield the 
DFE coefficients. Resultantly, the operation according to the present invention produces DFE 
coefficients faster and with less computational requirements &an the prior art operations, thus 
providing significant advantages over prior techniques. 

According to the present invention, a channel impulse response A is first estimated based 
upon either a known training sequence or an unknown sequence. A solution to the DFE 
coefficient computation problem is then cast as a standard recursive least squares (RLS) problem. 
More specifically, the solution for Feed Forward Equalizer (FFE) coefficients (of the DFE) 
based upon the channel impulse response h is formulated as the Kahnan gain solution to the RLS 
problem. A fast recursive method for computing the Kahnan gain is then directly used to 
compute g^. 

In one embodiment of the present invention, a fast transversal filter (FTF) technique is 
employed to compute g^^, . In this embodiment, the complexity of a conventional FTF algorithm 
is reduced to one third of its original complexity by choosing the length of Feed Back Equalizer 
(FEE) coefficients Z)„^ (of the DFE) to force the FTF algorithm to use a lower triangular matrix. 
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This technique significantly reduces the needed recursions and computations, as well as avoiding 
finite precision problems of a conventional FTF solution implemented in hardware. 

Finally, with the FFE coefficients g^^, determined, the FBE coefficients b^^, are 
computed by convolving the FFE coefficients g^^, with the channel impulse response k In 
performing this operation, a convolution matrix that characterizes the channel impulse response h 
is extended to a bigger circulant matrix. With the extended circulant matrix structure, the 
convolution of the FFE coefficients g^^, with the channel impulse response h may be performed 
the convolution operations in a transformed domain. In one embodiment, the convolution is 
performed in the firequency domain, which can be computed efficiently using the Fast Fouria: 
Transferal (FFT). However, in other embodiments, other transformation domains are employed 
for the convolution operations, e.g., Discrete Cosine Transformation domain and the Discrete 
Hadatnard Transform domain, among others. 

The method of the present invention is exact and is much faster than any prior technique 
that may have been employed to compute DFE coefficients for the same channel and DFE filter 
lengths. Because the method of the present invention relies upon simple structured recursion, the 
computation period and/or the hardware required for the DFE coefficient computations is 
significantiy reduced as compared to prior art implementations. The reduction in computation 
complexity and resultant increase in computation speed enables a DFE that operates according to 
the present invention to track fast channel variations. 

The method of the preset invention may be efficientiy used in many communication 
systems that require the use of equalizers such as mobile wireless receivers, fixed wireless 
receivers, cable modem receivers, HDTV receivers, etc. to increase the performance of such 
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communication systems. Other features and advantages of the present invention will become 
apparent from the following detailed description of the invention made with reference to the 
accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

These and other features, aspects and advantages of the present inv^tion will be more 
fully xmderstood when considered with respect to the following detailed description, appended 
claims and accompanying drawings wherein: 

FIG. 1 is a block diagram illustrating a discrete time symbol-spaced Decision Feedback 
Equalizer (DFE) channel equalization model 100; 

FIG. 2 is a logic diagram generally illustrating operation according to the present 
invention in determining DFE coefficirats and applying such coefficients to a DFE; 

FIG. 3 is a logic diagram illustrating operations according to the present invention 
employed to determine Feed Forward Equalizer (FFE) coefficients for the DFE; 

FIG. 4 is a logic diagram illustrating operations according to the present invention 
employed to determine Feed Back Equalizer (FBE) coefficients for the DFE; 

FIG. 5 is a block diagram illustrating a discrete time fractionally-spaced DFE that 
operates according to the present invention; 

FIG. 6 is a block diagram illustrating a multi-channel equivalent of the discrete time 
fractionally-spaced DFE of FIG. 5; and 

FIG. 7 is a block diagram illustrating a transceiver constructed according to the present 
invention. 
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DETAILED DESCRIPTION OF THE INVENTION 

FIG. 2 is a logic diagram generally illustrating operation according to the present 

invention in detennining Decision Feedback Equalizer (DFE) coefficients and in applying such 

coefficients to a DFE. The operations of the present invention are performed by a processor, 

such as a Digital Signal Processor (DSP), or o&er circuitry present within a receiver that 

determines DFE coefficients to be applied to a DFE, also resident in the receiver. The DFE 

operates upon samples of a received signal in an attempt to remove channel effects from the 

samples so that digital data may be extracted from the samples. The structure and operation of 

DFEs, one of which was illusfrated in FIG. 1, are generally known and will not be fiirther 

described herein except as they relate to the present invention. 

The operations of the present invention will be first performed by a processor, tap 
computer, DSP, or other receiver device, to determine initial DFE coefficients to be used in 
subsequent operations by flie receiver. Thus, during startup or reset, a channel corresponding to 
the receiver is estimated (step 202). According to one embodiment of the present invention, the 
channel is estimated based upon a known preamble sequence. However, in other embodiments, 
the channel could also be estimated based upon xmknown received data. In either case, channel 
estimation operations are generally well known and are not described further herein except as it 
relates to the present invention. 

With the channel estimated, Feed Forward Equalizer (FFE) coeffideats are determined 
based upon the channel estimate (step 206). Then, Feed Back Equalizer (FBE) coefficients are 
determined based upon the FFE coefficients and tiie channel estimate (step 208). The manner m 
which the FFE and FBE coefficients are generated is step 206 and step 208 will be described in 
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detail herein with reference to FIGs. 3-7. 

With the FFE and FBE coefficients determined, they are applied to the DFE (step 208) 
and are used in equalizing samples of a received signal to remove channel effects. These DFE 
coefficients are continually updated (step 210) using a known technique. Periodically, upon the 
receipt of a next packet for example, upon an indication that redetermination is required, or upon 
another triggering event, the DFE coefficients are redetermined (step 212). hi this event, another 
channel estimate may be obtained (step 214). Then, the DFE coefficients are again determined 
according to the present invention and applied to the DFE. The operations of FIG. 2 continue 
until the receiver is turned off, placed in a sleep mode, or otherwise mactivated. 

FIG. 3 is a logic diagram illustrating operations according to the present invention 
employed to determine Feed Forward Equalizer (FFE) coefficients for the DFE. hi a first 
operation of FIG. 3, a DFE delay is selected (step 302). hi one embodiment, this delay is selected 
as the channel length, hi such case, the DFE delay corresponds to the length of the FFE. 

Next, the DFE solution is formulated into a least squares solution (step 304). By collecting 
the FFE coefficients g and the FBE coefficients b mto a single vector w, the minimization of the 
Mean Square Error may be written as: 

Equation (9) 

Now, denoting R„ the variance of the augmented input regression vector u, and cross 
variance R^i„_s) > the well-known solution to this smoothmg problem is given by Equation (10) 
as: 



minE 



xin-S)-\y„-x„] 



10 



Attorney Docket No. BP 1914 



where 



Equation (10) 



A 



Equation (11) 



and 





-S)' 








-S) 







Equation (12) 



Using the channel output model of Equation (4), and the fact tiiat x(n) is individually 
identically disttibuted (i.i.d.), the following closed fonn expressions for 
{Ry,R^,R„R^^„-S)>RMn-s)\ are determined: 



R,= Ey:y„=<rU+H*H 

0, 



R^=H\^X) = H' 



0. 



A 



R^,„.sy=H*ExAn-S) = H'[0^ 1 Oj=h' 
where F is a submatiix of H as set forth in Equation (5), 



Equation (13) 
Equation (14) 

Equation (15) 

Equation (16) 
Equation (17) 



H = 



H 



Equation (18) 



11 



Attorney Docket No. BP 1914 

i?iis defined as the (<J + l)xL submatrix ofH, consisting of the first (^+1) rows ofH. Note 
that for the case of colored noise, the matrix cr,^/ should be replaced by the autocorrelation 
matiix of the noise . Extending the derivation to the colored noise case is stiaightforward. 
Now, with the quantities defined above. Equation (10) becomes: 



-H I 
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_0_ 



Equation (19) 



Using the well known inverse formula of block matrices, may be rewritten as: 
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and as: 
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{all^HlH.+HlH^r'h* 



which may be written as: 
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Equation (20) 



Equation (21) 



Equation (22) 



Althou^ each of the two matrices and has a shift structure, tiie augmented matrix 
does not have a shift structure. 

In selecting tiie length of the FBE (DFE delay) to beM ^ iV - 1 , the quantities hg , and 



are such that: 
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Hg 0' 
0 H 



Equation (23) 



and 



This implies that: 



h = {hs 0] 



Equation (24) 



<71i+h;h,+h;h,= 



all + HsH, 



0 



Equation (25) 



0 all+H'H 

is block diagonal. In this case, the expressions for and decouple into a simple form. 
Therefore, the optimal FFE and FBE coefficients are represented as: 



~s 



opt 



Equation (26) 
Equation (27) 



The above expressions are valid for all values of tibe DFE delay S . hi general, the optimal 
value for the DFE delay <y is within the range i-l<<J„^ <iV+X- 2. In the special case of the 

choice ^ = i - 1 , the matrices involved in the above expressions are given by 
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Equation (28) 
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Equation (29) 
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Note that this solution is equivalent to the solution of Al-Dhahir's equations, due to the 

uniqueness of when Equation (2) is minimized. However, the expressions obtained above 
for Equations (26) and (27) provide alternative methods to compute the FFE and FBE 
coefficients and in a simpler and more efficient way than the prior art techniques. 
In calculating g^^,, a coefficient matrix is first defined in Equation (30) as: 

= ((7^1 + H*Hs)-' Equation (30) 

so that the optimal solution for the FFE coefficients is given by 

g^=P,h;. Equation (31) 

The expression for g„j^ corresponds exactly to the definition of Hoe Kalman gain vector, 
used to update tiie optimal weights in a certain regularized Recursive Least Squares (RLS) 
problem. More specifically, given a (« + l)xZ data matrix H„ and the corresponding 
coefficient matrix P„ , the Kalman gain g„ = P„A* can be time-updated according to the following 
recursions: 

y-' (n) = 1 + h„P„_^h*„ , Equation (32) 

gn=Pn-Xm> Equation (33) 

P„ = Pn-x -gJ-\n)g: , Equation (34) 

where P_, = afl (the term P-i is the initial condition for P,0 and go = 0 • The computation of the 
Kalman gain vector g„^j in the above solution relies on the propagation of the Riccati 
variable P„ . This method of computation requires 0(1? ) operations per iteration. 
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Well known fast RLS schemes avoid the propagation of P„ and compute the gain ^„ in a 

more efificient way. In this case, the computational complexity required is only of 0{L) per 
iteration. This implies that the overall complexity needed to calculate the FFE coefficients can 
be on the order of O(i^) thus yielding an efficient method for computing the FFE coeffideots. 
Therefore, fast RLS filters are used to determine . hi such case, fast transversal computations 
are employed to determine the Kahnan gam for the RLS solution (step 306). Here we also note 
tfiat it is straightforward to extend tfiis method to use other fast RLS algorithms, e.g., array form 
fast RLS algorithms. 

Faster recursions may be obtained by selecting the FBE length. Such selection eliminates 
certain variables that remain constant and equal to zero during adaptations in the special solution 
for g^. Fast RLS recursion in its explicit form propagates g„ in an efficient manner, see, e.g.: 
[1] Ljung, M. Morf, and D. Falconer, "Fast calculation of gain matrices for recursive estimation 
schemes," hit J. Contr. vol. 27, pp. 1-19, Jan 1978); [2] G. Carayannis, D. Manolakis, and N. 
Kalouptsidis, "A fast sequential algorithm for least squares filtering and prediction," IEEE Trans, 
on Acoustic, Speech, Signal Proc, vol. ASSP-31, pp. 1394-1402, December 1983; and [3] J. 
CiofB and T. Railath, "Fast recursive-least-squares transversal filters for adaptive filtering," 
IEEE Trans, on Acoust., Speech Signal Processing, vol. ASSP-32, pp. 304-337, April 1984. 

Table 1 lists the fast recursions employed in one embodiment for computing the normalized 
Kahnan gain. The additional index included in kj^^ and fj^in) stems from the fact that these 
quantities admit an order-update relation, instead of time-update. 
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Initialization 

1 = = = 0 

n(o)=i 

Forn = 0to S repeat operations fl) to (13): 

(1) a^{n-\) = h{n)-h„_y„_^ 

(2) f{n-\) = n{n-\)a{n-\) 



0 


_ a*{n-\) 


1 


_^I,n-l _ 







(4) C^{n-\) = C^{n-l) + a*{n- \)f{n - 1) 

(5) >^-, = ^-2+^i.«-i/(«-l) 

(6) /..i(/^) = n("-l) '^/^^_^^ 

(7) v(n) = (last entry of Ar^, „_i ) 

(8) = ^L^in-i 0 : + v(«)w*_i 

(9) A«) = r(«-l)v(») 

(11) = 

(12) C'in) = C'{n-l) + fi\n)bin) 

(13) w:^=w*_i+^,X«) 

Table 1: Fast Transversal Computation of the Kalman gain. 

The purpose of computing k„ in the well-known FTF algorithm is to use it in the 
computation of the corresponding optimal least squared solution for the FFE coefficients. 
Because we are interested only in k„ , the filtering part of the FTF algorithm is not necessary, and 
does not appear in the algorithm listing. 

The quantities are known as the least-squares solutions of Ihe forward and 
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backward prediction problons of order i, with corresponding residual errors {fin),b(n)}- Now, 
because during the first L iterations, the desired signal for the backward prediction problem is 
equal to zero, the backward least squares solution vi^ will be equal to zero. This means that all 
quantities associated with the backward prediction problems will remain nulls for the first L 
iterations. Since in our case the optimal solution is achieved exactly at the L-th iteration, we can 
simply rule out the computation of these quantities Ifrom Table 1. These operations correspond to 
operations (7) and (9) through (13) of Table 1. Note that operation (10) of Table 1 is also 
eliminated, since /3(n) = 0, which implies yi(n) = ri+i(«)' This i^^eans that we can replace 
operation (6) of Table 1 with: 



which can be further simplified since 



Moreover, operation (8) of Table 1 becomes simply 



Equation (35) 



= m EquationOe) 
^ CUO) CUV) CUn-l) CUn-l) 



Equation (37) 



where (1:L) denotes the first L entries of A;i+i,„_i 

Table 2 illustrates a simplified fast recursion for computing tiie optimal FFE coefficients 

S opt * 



17 



Attorney Docket No* BP 1914 



Initialization 



r(-i)=^ 

KO)=i 



For n = Oto d repeat operations (1) throu2h (7): 
(l)a,{n-\) = h{n)-K_y„_, 
0) f{n-\) = y{n-\)a{n-\) 



(3) Wi = 



0 


, /*(«-!) 


1 









(4) ^^(« - 1) = ^^(« - 2) +«* (« - !)/(« - 1) 

(5) h/,=m/2+^,,„_/(«-1) 

(6) r(«) = 



Table 2: Fast Transversal Computatioii of the FFE coefficients. 
With the Kalman gain determined, the FFE coefficients are tiien determined (step 308). In 
tiie recursions of Table 2, the FFE coefficients are determined by setting = kj^^yiS) when the 

number of iterations, n, is equal to the DFE delay. 

Unlike the conventional weighted FTF algorithm, the above recursions do not face finite 
precision difficulties for the following reasons. First, by ruling out the equations associated with 
the backward prediction problem, we are automaticdly eliminating many of the recursive loops 
that are responsible for the finite precision difficulties of the fiill FTF algorithm. Second, tiiese 
simplified fast recursions have a forgetting factor A = l, which yields finite precision stability. 
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Third, the simplified algorithm deals with a finite set of data (S + 1 iterations). This algorithm is 
tiien reset, which avoids the accumulation of finite precision errors. 

FIG. 4 is a logic diagram illusteating operations according to the present invention 
employed to determine Feed Back Equalizer (FBE) coefficients for the DFE. The FEE 
coefficients b^p^ are determined according to the matrix-vector product of Equation (27). The 
computation of the feedback filter coefficients simply amounts to convolution of the channel 
impulse response with g^p^ . The convolution operation that defines the optimal FBE coefficients 
in Equation (27) can be computed directiy with LM/2 multiplications. Alternatively, the 
operations may also be efficiently performed using well-known fast FFT convolution techniques. 
To illusti-ate such operations. Equation (29) is first rewritten in Equation (38) as: 



'h 


Hi 




Sopt 








0 



Equation (38) 



The above equality holds true regardless of the values of H^, H^, or . Now, if these values 
are chosen such that the matrix C, defined by: 



C = 



H 



Equation (39) 



is a square (Q x Q) circulant matrix, where Q is the smallest power-of-two integer larger than or 
equal M+ In this case, the matrix C is rewritten in Equation (40) as: 

C = F * AF , Equation (40) 

where F is a (gx0 FFT matrix and A is a diagonal matrix tfiat contains the elements of the 

FFT of the first row of C. The solution for b^^^ becomes: 
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Sopt 

0 



Equation (41) 



The complexity of this method is the complexity of obtaining the FFT of two vectors of Q 
5 elements of each, the inverse FFT of another Q element vector, and Q complex multiples. Thus 
the overall complexity is 

Q+3Qlog2iQ) ■ Equation (42) 

1 0 For the case of a power-of-two channel estimate N, the complexity is 

2M + 6Mlog2(2M). Equation (43) 

b Thus, referring again to FIG. 4, in determining the FBE coefficients. Hie convolutional 

Q _ 
1^ matrix His first determined (step 402). As was previously described, the matrix H may be 

Z determined as part of the channel estimation. Then, the convohitional matrix H is extended to 

W 

^ form the bigger circulant matrix C (step 404). The biggo" drculant matrix C is then converted to 
P the Frequency domain (step 406) as are the FFE coefficients (step 408). With both the bigger 
2 drculant matrix C and the FFE coefficients in the firequency domain, they are convolved in the 
20* frequency domain by simple matrix multiplication (step 410) to produce the FBE coefficients 
b„^. Finally, the resultant FBE coefficients are then converted to the time domain using 
inverse FFT operations to produce FBE coefficients b^p, in the time domain (step 412). The FBE 
coefficients are then q)plied with the FFE coefficients to the DFE (as previously 



described with reference to FIG. 2. 
25 FIG. 5 is a block diagram illustrating a discrete time firactionally-spaced DFE that 

operates according to the present invention. FIG. 6 is a block diagram illustrating a multi- 
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channel equivalent of the discrete time fractionally-spaced DFE of FIG. 5. The approach xised 
for the symbol-spaced model can be easily extended to fractionally spaced models, hi this 
section, we shall derive a fast algorithm for a T/i-spaced equalizer. The fast solution for flie 
general case of faster oversampling factors, say T/i, T/4, . . .> etc. will then follow simply by 
mspection of the arguments in this section. 

FIGs. 5 and 6 illustrate two equivalent representations for a discrete time model of a 272- 
spaced equalizer. FIG. 5 illustrates a DFE 500 in terms of an down sampler and up sampler 
device, while FIG. 6 illustrates a DFE 600 is a corresponding multichannel representation. Thus, 
as contrasted to FIG. 1, the DFE 500 of FIG. 5 includes a down sampler 502 and an up sampler 
504 that alter the sampling frequency. Further, the DFE 600 of FIG. 6 includes a first channel 
function 602 having channel estimate \ and corresponding noise fimction Vq(«) and FFE 606 

coefficients . Likewise, the DFE 600 of FIG. 6 also includes a second channel fimction 604 

having channel estimate \ and corresponding noise fimction {n) and FFE 608 coefficients 

The equivalence between the structures of FIG. 5 and FIG. 6 can be easily verified by 
writing each of the quantities {^.g} and v(«) in terms of their polyphase components, and 
interchanging their position with the down sampler and the up sampler (see, e.g., P.P. 
Vaidyanathan, "Multirate Systems and Filter Banks," Prentice Hall, NJ, 1993 and J.R. Treichler, 
1. Fijalkow, C.R. Johnson, Jr., "Fractionally spaced equalizers," IEEE Signal Processing 
Magazine, vol. 13, no. 3, May 1996 for details). The quantities {Ao,/ri} and {go?5i} ^re the so- 
called polyphase components of the channel and FFEs, and are denoted by N and L size vectors 
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{h^.h^} and {go.gi} (see FIG. 6). The noise sequences {vo(n\v^(n)} are also the even and odd 
samples of v(n), with powers 2(T^ . 

A model for the multichannel representation, just like the one for the symbol spaced case in 
equation (4) may be rewritten by collecting {go?^i} i^to a single vector g' , as: 



/=[^i go\ 

with the output of the multichannel system given by: 

Thus, the following model of Equation (46) works for the input to g' : 



Equation (44) 



Equation (45) 



k« yunl^^nlHo ^i]+k„ Equation (46) 

where {Hq,H^} are the convolution matrices associated with the subchannels {/?o>^}- K is more 

convenient, however, to express the inner product in Equation (45) as a function of tiie original 
FFE coefficients vector g, by reordering the entries of . In this case, Equation (46) is replaced 
by: 



yn=^nH +V„ 



where 



h(3) 



A(0) 
h(2) 



h(3) 



HO) 

h(2) 



h(2N-l) h(2N-2) 



h(2N-\) h(2N-2) 



Equation (47) 



KO) 
h(2) 

h{2N-\) h(2N-2) 



Equation (48) 
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Given this convolution matrix and tiie noise variance <j^. , the solution for the fractionally space 
problem is simply given by Equations (26) and (27), with {fir^,F,cr^} replaced by \Hs,H\al\ 
(similarly to the symbol-spaced case, here we assume M >N-1). 

For the symbol-spaced DFE, the shift structure of implies that A:^ „ = ^i:+i,„-i (1 ^) • 
That is, the normalized gain k^ ^ can be computed fast, by performing the order update of k^ „_^ 
to , and then retaining the first L entries of as A:^ „ . The procedure to derive a fast 

recursion in the r/2-spaced case follows Hie same principle* The only difference here is that the 
data matrix H'^ has now a double shift structure, in the sense that each row is formed by shifting 

two samples of the channel at a time. Thus, two successive order updates from k^ „_^ to A:^+2,«-<i 

are performed and then the first L entries of A:^+2^-i retained: 

ki^n-i ^x.2,«-i Equation (49) 

so that 

= ^i+2,«-i (1 • Equation (50) 

In other words, the resulting algorithm is given by two consecutive forward prediction problems 
similar to operations (1) through (6) of Table 2 of orders L and Z + 1, respectively. Table 3 lists 
the resulting algorithm for the 772-spaced case. 
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InitiaUzation 




= 0 



For rf^ Oto S repeat operations (I) to (II): 



(1\ For q^Lto L+1 repeat operations (I) to (6): 



(1) (n - 1) = h(2n + q-L + l)- h„_y,^_2 



0 




-1) 


1 






-2) 





(4) C(«-l) = C("-2)+<(«-l)/,(«-l) 

(5) <„-i =wi„-2 +^,,„-i/,(«-l) 



Table 3: Fast Transversal Computation of FFE coefficients for the r/2-spaced equalizer. 
In the general case of a 775 spaced equalizer, the fast DFE tap computation is a 
straightforward extension of the above algorithm, where q is iterated from L to (L+S-l). 
However, in the fractionally spaced equalizer, the RLS problem is formulated as a multi-channel 
problem. In such case, a multi-channel Kalman gain is calculated for the multi-channel RLS 
problem and the FFE taps are determined therefrom. Note that in the expressions of Table 3, 
successive order updates are performed. 



(6) r,+i («) = /,(«) 



(11) 



(l:L) 



Set g,p,=k, srLAS) 
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Now, the optimal FBE taps could now be computed as 



KptAlM 0]F*A,(F®/,) 




Equation (51) 



where is a block diagonal matrix, which satisfies 



Equation (52) 



where Fs-F®!^ and is a block circulant matrix, which is formed by extending the matrix 
if in the same manner as in Equation (39). As was the case with the symbol spaced operations 
described above, the convolution may be performed in a transformed domain, e.g., Frequency 
Transformation domain. Discrete Cosine Transformation domain and the Discrete Hadamard 
Transformation domain, among others. In such case a Croniker product may be employed in 
conjunction with such domain transformation in performing multi-channel convolution in the 
selected transformed domain. 

FIG. 7 is a block diagram illustrating a transceiver constructed according to the present 
invention. The components of the transceiver 700 are resident in a communications device and 
are illustrated only generally to show how liie operations of the present invention would be 
accomplished in such a transceiver. The transceiver 700 includes a receiver section 702 and a 
transmitter section 704 and further includes either a wireless transceiver 706 or a wired 
transceiver, depending upon the system in which the transceiver 700 is implemented. In the case 
of a cellular device, RF device, satellite system device, or another wireless implementation, the 
transceiver 700 includes a wireless transceiver. However, in the case of a cable modem, a LAN 
device, a home networking device, or another device that couples to a physical media, the 
transceiver 700 includes a wired transceiver 708. Further, if the present invention is 
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implemented in a Serializer/Deserializer (SERDES) or similar application, the receiver section 
702 and transmitter section 704 may couple to a media without a wired transceiver 708. 

Further discussion of the transmitter section 704 and receiver section 702 are in the 
context of baseband processing. In such case, the receiver section 702 receives a baseband signal 
from the wireless transceiver 706 (or wired transceiver 708) that is baseband modulated and 
operates upon the baseband signal to extract data. These operations include determining DFE 
coefficients according to the present invention and operating upon the baseband signal using the 
determined DFE coefficients. 

The transmitter section 704 receives digital data to be transmitted from a host, codes the 
digital data into a baseband signal, and passes the baseband signal to the RF transceiver 706. The 
RF transceiver 706 couples the baseband signal to an RF carrier to create an RF signal and 
transmits the RF signal to a receiving device across a wireless link. 

The receiver section 702 receives a baseband signal that carries coded data from the RF 
transceiver 706. A Programmable Gain Amplifier (PGA) 712 adjusts the gain of the baseband 
signal and then provides the gain-adjusted baseband signal to an Analog-to-Digital Converter 
(ADC) 714 for sampling. The ADC 208 samples the gain adjusted baseband signal at a 
particular sampling frequency, fg (that is the symbol clock frequency), to produce samples fliereof. 

A processor 710 couples to the output of the ADC 714 and analyzes a preamble sequence 
contained in each received physical layer frame. Based upon the preamble sequence, the 
processor 710 determines a gain to be applied to portions of the baseband signal corresponding to 
the data carrying portions of the frame and provides this gain to the PGA 712. Further, the 
processor 710 may also interact with the optional timing compensation section 716 to 
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compensate for symbol timing and RF carriCT mismatches. 

The processor 710, based upon the preamble sequence (and based upon actual extracted 
data in some operations), also determines FFE 104 and FBE 108 coefficients. The manner in 
which these coefficients are determined was previously described in detail herein. Further, and 
as was also previously described, the processor 710 may estimate a channel and calculate DFE 
coefficients based upon unknown but assumed data content. After the processor 710 determines 
these coefficients, they are applied to the FFE 104 and FBE 108 for subsequent use in extracting 
data from the baseband signal. 

The structure described in FIG. 7 may be embodied using various types of circuits formed 
using various manufacturing processes. For example, in one particular embodiment, the 
transceiver 706 (or wired transceiver 708) is embodied in a first integrated circuit that is coupled 
to a second integrated circuit that includes the transmitter section 704 and the receiver section 
702, among other circuits. In another embodiment, the RF transceiver 706, the transmitter 
section 704 and the receiver section 702 are all formed on a single monolithic integrated circuit. 
These integrated drcuits may be constructed in CMOS or another semiconductor technology, 
e.g., PMOS, NMOS, Bipolar, etc. 

Further, the receiver section 702 of FIG. 7 may be constructed using various circuit 
elements/combinations, hi one embodiment, all structures past the ADC 714 m the receiver 
section 702 are embodied using a Digital Signal Processor (DSP) or similar processing device. 
In another embodiment, dedicated signal path circuitry embodies each of the structural 
components of the receiver section 702, including the processor 710. While a DSP 
implementation would provide more flexibility, dedicated signal path circuitry would typically 
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provide hi^er performance at a lower cost and with lower power consraaption. 

The structure and operation of the present invention may be employed in satellite 
communication systems, satellite television systems, HDTV systems, fixed wireless 
communication systems, mobile wireless communication systems, cable modem/television 
communication systems, home networking systems, wireless local area networking systems, 
wired local area networking systems, and many other types of communication systems. The 
present invention applies to all types of communication devices in which equalizers are 
employed to operate upon received signals. 

The invention disclosed herein is susceptible to various modifications and alternative 
forms. Specific embodiments therefore have been shown by way of example in the drawings and 
detailed description. It should be understood, however, that the drawings and detailed 
description thereto are not intended to limit the invention to the particular form disclosed, but on 
the contrary, the invention is to cover all modifications, equivalents and alternatives felling 
within the spirit and scope of the present invention as defined by the claims. 
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